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We demonstrate that the Configuration Interaction (CI) Approximation recaptures essential 
features of the exact (Bethe Ansatz) solution to the ID Hubbard model. As such, it provides valu- 
able route for describing effects which go beyond mean-field theory for strongly correlated electron 
systems in higher dimensions. The CI method systematically describes fluctuation and quantum tun- 
neling corrections to the Hartree-Fock Approximation (HFA) . HFA predicts that doping a half-filled 
Hubbard chain leads to the appearance of charged spin-polarons or charged domain-wall solitons 
in the antiferromagnetic (AFM) background. The CI method, on the other hand, describes the 
quantum dynamics of these charged magnetic solitons and quantum tunneling effects between var- 
ious mean-field configurations. In this paper, we test the accuracy of the CI method against the 
exact solution of the one-dimensional Hubbard model. We find remarkable agreement between the 
energy of the mobile charged bosonic domain-wall (as given by the CI method) and the exact energy 
of the doping hole (as given by the Bethe Ansatz) for the entire U/t range. The CI method also 
leads to a clear demonstration of the spin-charge separation in ID. Addition of one doping hole to 
the half-filled antiferromagnetic chain results in the appearance of two different carriers: a charged 
bosonic domain- wall (which carries the charge but no spin) and a neutral spin-1/2 domain wall 
(which carries the spin but no charge). 
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I. INTRODUCTION 

The two-dimensional (2D) Hubbard model (or a modi- 
fied version of it) is widely believed to describe the essen- 
tial physics of the strongly correlated electrons in high 
temperature superconducting cuprates.a Unlike the ID 
Hubbard model, an exact solution of the 2D problem 
is not known, and one must resort to approximations. 
The validity of such approximations can be tested on the 
ID Hubbard model, which is exactly solvablea using the 
Bethe Ansatz E It is highly desirable to develop an ap- 
proximation which recaptures essential physical features 
of the exact Bethe Ansatz solution and which, at the 
same time, can be applied to higher dimensional systems. 
In this paper-, we show that the Configuration Interaction 
(CI) Methodel is such an approximation. We demonstrate 
that the CI method recaptures essential quantum tunnel- 
ing effects that go beyond mean-field theory and that lead 
to spin-charge separation in ID. The predictions of the CI 
method in the 2D case as well as a comparison to cuprate 
superconductors are described in detail elsewhereB 

We consider the generalized Hubbard Hamiltonian: 

n = -* E { c lT%c j0 + h.c) + uj2 4^4i c *i (!) 

(i.j) i 
a/3 

Here, cj is the creation operator for an electron of spin 
a at site i, and the notation means that the sum is 
restricted to nearest- neighbors sites i and j. The param- 
eters of the problem are the hopping matrix t, the on- 
site interaction matrix U, and the SU(2) matrices TMj, 
which describe phase factors and internal spin rotations 
acquired by the electron as it hops between sites. For the 
ID Hubbard Hamiltonian which we consider in this paper 
we simply set T % K — S a/ 3. However, for the 2D system, 
it has been shown that a non-trivial choice of the matri- 
ces Tja, which describes a 2ir (internal) spin rotation of 
the electron as it encircles an elementary plaquette of the. 
square lattice, may be essential to the cuprate physicsfm 
For this reason, we develop the CI formalism for the more 
general Hamiltonian ([[]). 

In addition to providing a test to gauge the accuracy 
of the Configuration Interaction Method, the ID results 
may be directly relevant to certain high T c cuprate su- 
perconductors. YBa2Cu307 and its close relatives have 
quasi-one-dimensional CuO chain structures. Experi- 
ments measuring the dc resistivityQ, the infrared and 
optical conductivity^ and the pepetration depth in un- 
twinned crystalal3 and ceramics^ have revealed large 
anisotropics between the a-direction (perpendicular to 
chains) and the b-direction (parallel to chains). These 
results suggest that substantial currents are carried along 
the chains in both the normal and superconducting state. 
The source of superconducting condensate on the chains 
has not yet been elucidated. 

We begin Section 2 with a brief review of the Static 
Hartree-Fock Approximation (HFA). The HFA leads to a 



mean-field ground-state with properties which are in dis- 
agreement with those of the exact ground-state. For the 
half-filled (undoped) chain, HFA predicts a degenerate 
ground-state with long-range antiferromagnetic (AFM) 
order and with staggered magnetic moments aligned 
along some arbitrary direction. Even in the presence of a 
weak external interaction which creates an easy axis for 
spin orientation, the HF ground state is degenerate with 
the related mean-field in which all the spins have been 
flipped (S — > — S). On the other hand, the true ground- 
state of the half-filled system is non-degenerate and has 
no long-range order, despite the presence of strong AFM 
correlations. When holes are introduced into the chain, 
the HFA leads to the creation of static magnetic spin- 
polaron and domain-wall solitons which trap the hole. 
Since these solitons are motionless in the static HFA, this 
mean field solution breaks the translational invariance of 
the original Hamiltonian. The CI method facilitates the 
restoration of both of the above symmetries in accord 
with the exact Bethe ansatz ground state. 

The essence of the CI method is to use a linear combi- 
nation of HF wavefunctions in order to restore the various 
broken symmetries of the mean-field theory. In a doped 
system, for instance, the CI wavefunction is chosen to be 
a linear combination of HF wavefunctions describing the 
charged soliton centered at various sites. Besides restor- 
ing the translational symmetry, such a wavefunction also 
effectively takes into account the quantum dynamics of 
the charged soliton along the chain. This motion repre- 
sents a large amplitude tunneling event between a given 
AFM mean field and the alternative AFM mean-field ob- 
tained by the operation (S — > —5). Moreover, the prop- 
agating soliton lowers its quantum zero point energy con- 
siderably, relative to the static (HFA) soliton. 

In Section 3 we briefly review the exact Bethe Ansatz 
(BA) solution and the ground-state energies of the half- 
filled and doped chains. The comparison between the 
Bethe Ansatz, Hartree-Fock and Configuration Interac- 
tion solutions is presented in Section 4. For both the 
undoped chain, and the chain with one doping hole, we 
show that through the nucleation of mobile quantum soli- 
tons, the CI method provides a much better description 
than the HFA. We identify the charge carrier of the Hub- 
bard chain to be a charged bosonic domain-wall, and we 
show that its energy, in the CI method, is in excellent 
agreement with the BA result. The CI method leads to a 
simple, physical interpretation of spin-charge separation. 
Adding a single hole to a chain leads to the creation of 
two magnetic solitons. One soliton is the bosonic charged 
domain wall (which carries charge, but no spin), while the 
other soliton is the neutral spin- 1/2 domain wall (which 
carries spin, but no charge). The energy of these excita- 
tions are in good agreement with the BA results. Finally, 
Section 5 contains discussion of the results and conclu- 
sions. 
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II. APPROXIMATIONS: HARTREE-FOCK AND 
CONFIGURATION INTERACTION METHOD 

A. The Static Hartree-Fock Approximation 

One of the most widely used approximations for the 
many-electron problem is the Static Hartree-Fock Ap- 
proximation (HFA). In this approximation the many- 
body problem is reduced to one-electron problems in 
which each electron moves in a self-consistent manner 
depending on the mean-field potential of the other elec- 
trons in the system. While this method is insufficient, 
by itself, to capture all of the physics of low dimensional 
electronic systems with strong correlations, it provides a 
valuable starting point from which essential fluctuation 
corrections can be included. In particular, we use the 
Hartree-Fock method to establish the electronic struc- 
ture and the static energies of various magnetic soliton 
structures. In the more general Configuration Interaction 
(CI) variational wavefunction, the solitons acquire quan- 
tum dynamics and describe large amplitude tunneling 
and fluctuation effects that go beyond mean field theory. 

In the HF approximation, the many-body wavefunc- 
tion <3>) is decomposed into a Slater determinant of ef- 
fective one-electron orbitals. The one-electron orbitals 
are found from the condition that the total energy of the 
system is minimized 



5 W =0 



(2) 



In order to approximate the ground state of the 
Hamiltonian (I), we consider a Slater determinant trial- 
wavefunction of the form 



i*>=rRi°>' 



(3) 



where |0) is the vacuum state, N e is the total number of 
electrons in the system and the one-electron states are 
given by 



4 = 5^0"(*,cr)ct CT 



(4) 



Here, the one-particle wave- functions <p n (i,a) form a 
complete and orthonormal system. 

Using the wavefunction (j^) in equation (||), and min- 
imizing with respect to the one-particle wavefunctions 
cj) n (i,cr), we obtain the Hartree-Fock eigen-equations: 

E n <j> n (i,a) = -t ^2 T «/3^n(i,/3) 
+ uY / (b a pQ(i)-v*pS(i))Mi,[3) (5) 



where (cr x , u y ,a z ) are the Pauli spin matrices and the 
charge density, 



qwh*i4*«i*> = £i#p(<>«)i s 



(6) 



and the spin density, 



must be computed self-consistently. The notation j £ 
Vi appearing in (||) means that the sum is performed 
over the sites j which are nearest-neighbors of the site i. 
The self-consistent Hartree-Fock equations must 
be satisfied by the occupied orbitals p = 1 . . . N e , but can 
also be used to compute the empty (hole) orbitals. 

The ground-state energy of the system in the HFA is 
given by 



E GS = {mm =J2 E v~ u H ( i^w 2 - w 

p—l i ^ 



(8) 



where the single particle energies are obtained from (|j). 

The approximation scheme described so far is called 
the Unrestricted Hartree-Fock Approximation, because 
we did not impose constraints on the wavefunction 
which would require it to be an eigenfunction of various 
symmetry operations which commute with the Hamilto- 
nian ([I]) . If these symmetries are enforced, the method is 
called the Restricted Hartree-Fock Approximation. We 
use the Unrestricted HFA since it leads to lower ener- 
gies. The breaking of symmetries in our case implies 
that electronic, correlations are more effectively taken 
into account .113 The restoration of these symmetries is 
deferred until the CI wavefunction is introduced. 

In the undoped (half-filled) case, the self-consistent 
Hartree-Fock equations can be solved analytically for 
the infinite system, using plane-wave one-particle wave- 
functions. In the unrestricted Hartree-Fock approach, 
doping the system leads to the appearance of inhomoge- 
neous solutions, which break the translational invariance. 
In this case, we solve the unrestricted self-consistent 
Hartree-Fock equations numerically on a finite chain. 
Starting with an initial spin and charge distribution S(i) 
and Q(i), we numerically solve the eigenproblem (^) and 
find the HF eigenenergies E n and wavefunctions 4> n (i, a). 
These are used in Eqs. (||) and (Q) to calculate the new 
spin and charge distribution, and the procedure is re- 
peated until self-consistency is reached. Numerically, we 
define self-consistency by the condition that the largest 
variation of any of the charge or spin components on any 
of the sites of the lattice is less that I0~ 9 between suc- 
cessive iterations. 
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B. Configuration Interaction Method 

The basic idea of the Configuration Interaction (CI) 
method is that the ground-state wavefunction, for a sys- 
tem with N e electrons, is not just a N e X N e Slater deter- 
minant (as in the HFA), but a judiciously chosen linear 
combination of such Slater determinants. □ Given the fact 
that the set of all possible Slater determinants (with all 
possible occupation numbers) generated from a complete 
set of one-electron orbitals constitute a complete basis of 
the iV e -particle Hilbert space, our aim is to pick out a 
subset of Slater determinants which captures the essen- 
tial physics of the exact solution. 

Consider the CI ground-state wavefunction given by 



N 



|*)=^a i |* i ) 



(9) 



where each is a distinct N e x N e Slater determinant 
and the coefficients 014 are chosen to satisfy the minimiza- 
tion principle: 



_6_ f(V\H\V) 

da. 



= 



(*l*> 

This leads to the system of CI equations 

N 



i = l,N 



(10) 



N 



^ Hijdj = E ^ °ij a j i = l,N (11) 



where E = (^\TL\^i) / (^\^) is the energy of the system in 
the \*f?) state , Ttij = (^ilH^j) are the matrix elements 
of the Hamiltonian in the basis of Slater determinants 
= l,N}, and = are the overlap ma- 

trix elements of the Slater determinants (which are not 
necessarily orthogonal). The CI solution is easily found 
by solving the linear system of equations ( |ll"l) , once the 
basis of Slater determinants — l,N} is chosen. 

If we denote by 0£ '(i,cr) the p — l,...,N e one-electron 
occupied orbitals of the Slater determinant \^ n ), these 
matrix elements are given by: 



O — 



frnra 



PI 



■n m 



P 



■ii in 

NeA 



P 



nrn 
N e ,N e 



(12) 



The matrix elements of the Hamiltonian (]]]) can be writ- 
ten as: 



(13) 



where the expectation values of the hopping and on-site 
interaction terms are: 



N 
p=l 



P 



II III 
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Here. 



E 

<i,3> 
a0 



,(<) = C ) *(<t)C ) (it) 



We now consider the specific choice of the Slater de- 
terminant basis {\^i),i — 1,N}. Strictly speaking, one 
may choose an optimized basis of Slater determinants 
from the general variational principle: 



(*l*> 







l,N;p=l,N e . (14) 



However, implementation of this full trial-function min- 
imization scheme (also knowiij-as a multi-reference self- 
consistent mean-field approachta) is numerically cumber- 
some even for medium-sized systems. Instead, we select 
the Slater determinant basis = 1,N} from the 

set of broken symmetry, Unrestricted Hartree-Fock wave- 
functions (||), their symmetry related partners and their 
excitations. Clearly (|^) satisfies (Q) by itself, provided 
that the a coefficients corresponding to the other Slater 
determinants in Eq.(|9|) are set to zero (see Eq. (0)). Since 
this Unrestricted HF wavefunction is not translationally 
invariant (the doping hole is always localized somewhere 
along the chain), we can restore the translational invari- 
ance of the CI ground-state wavefunction by also includ- 
ing in the basis of Slater determinants all the possible 
lattice translations of this Unrestricted HF wavefunction. 
In 2D, we must also include all the possible non-trivial 
rotations. 

Clearly, all the translated HF Slater determinants lead 
to the same HF ground-state energy (^ n |W|^ n ) = Eqs 
as defined by Eq. (||). The CI method lifts the degen- 
eracy between states with the hole-induced configuration 
localized at different sites, thereby restoring translational 
invariance. We may identify the lowering in the total 
energy due to the lifting of this degeneracy as quan- 
tum mechanical kinetic energy of deconfinement which 
the doping-induced configuration saves through hopping 
along the chain. In addition, quantum fluctuations in 



4 



the internal structure of a magnetic soliton can be incor- 
porated by including the lowest order excited state con- 
figurations of the static Hartree-Fock energy spectrum. 
Such wavefunctions are given by a^ah\^), where p > N e 
labels an excited particle state and h < N e labels the 
hole which is left behind (see Eq. (j^)). Once again, 
all possible translations of this "excited" configurations 
must be included in the full CI wavefunction. These ad- 
ditions can describe changes in the "shape" of the soliton 
as it undergoes quantum mechanical motion through the 
crystal. 



III. EXACT SOLUTION OF THE ID HUBBARD 
MODEL: THE BETHE ANSATZ 

Before reporting the results obtained in the HF and CI 
approximation for the ID Hubbard model, we briefly de=. 
scribe the exact Bethe Ansatz solution of this problem,cl 
for comparison purposes. Consider an A^-site chain with 
N e electrons of .which M have spin-down. Here, N e < N 
and M < N e /2a Any wavefunction satisfying the many- 
body Schrodinger equation 7i\^f) = E^) is character- 
ized by N e quasi- momenta — ^ < kj < ^ describing 
the motion of the N e electrons (a is the lattice con- 
stant), and M rapidities, A Q —describing the spin order- 
ing. Using the Bethe AnsatzH (BA) and imposing peri- 
odic boundary conditions, it can be shown that the quasi- 
momenta and th6..*apidities satisfy the so-called Bethe 
Ansatz equationsatj 



M 



exp (ikjNa) = J^J 



sin (kjd) — A a + Iji 
~ x sin (kjd) — A Q — ^-i 



(15) 



N = N e c 
be given b 



in the thermodynamic limit is known to 



E = -A\t\N 



dxJ (x)Ji (x) 



o x(l +exp(xU/(2|t|))) 



(18) 



where the J v (x) are cylindrical Bessel functions. In Fig. 
|l| we plot the ground-state energy per site, in units of t, 
obtained for a half-filled chain of various lengths N , for 
U/t = 5. While for very low values of N there are large 
variations between the energies of chains with even and 
odd number of unit cells, as N increases the energies ob- 
tained converge towards the thermodynamic value of Eq. 
( |ilf ) (shown as the full line). We consider chains with an 
even number of sites (integer number of unit cells), since 
we know that the ground-state has AFM correlations. 
Similar curves are obtained for other values of U/t. 

Another well known case is that of a chain with just 
one hole, in the U/t — > oo limit. In this limit, double oc- 
cupancy is forbidden by the large on-site interaction, and 
the only possible charge fluctuation is the motion of the 
hole. The total energy of the chain reduces to the energy 
of the hole, since the single occupied sites give little con- 
tribution to energy (in this limit the contributions from 
the AFM correlations of the electron spins, of the order 
of t 2 /U, are negligible for finite chains). It is straightfor- 
ward to showEd that the hole's dispersion relation in ID is 
exactly that of a free particle, and therefore the total en- 
ergy of the chain containing a single hole with momentum 
hk is E(k) — —2tcos(ka). In the following sections, we 
show that the ground-state energy approaches E = — 2t 
as U/t — > oo for the CI method as well as for the Bethe 
ansatz equations of the chain with one hole. 



Yj A Q — sin (kjd) + 



n 

0=1 



A„ - kr- 



A„ - A 







(16) 



The total energy and the total crystal momentum are 
then given by 



JV,. 



E = —2 y~] cos (fcj-a) 



P = J> (17) 



The ground state is always given by real fc's and A's. Ex- 
cited states are usually described by complex rapidities 
in so-called "string" structures. □ 

We solve the Bethe Ansatz equations iteratively, start- 
ing with a guess for the set of real rapidities A a (related 
to the ground-state solution of the ID Heisenberg chain, 
as described in ReferenceE3). Then, we solve Eq.( |15|) and 
find the quasi- momenta kj, which we use in Eq.(|16|) to 
find the new set of real rapidities. The procedure is re- 
peated until self-consistency is reached. 

We can check our numerical procedure in two partic- 
ular cases. First, the ground-state energy of a half-filled 



IV. COMPARISON BETWEEN THE BA, THE CI 
AND THE HF RESULTS 

A. The undoped ground-state 

1. Hartree-Fock results 

For the undoped system, the self-consistent HF equa- 
tions (^,^|,0) for an infinite system give rise to a mean- 
field ground state with long range antiferromagnetic or- 
der. However, the Mermin- Wagner theorem states that 
the true ground-state of any one dimensional isolated sys- 
tem cannot have long range order (LRO) and that LRO 
is absent in 2D systems for non-zero temperatures. In 
the framework of the CI method, mobile solitons in the 
AFM background mediate the destruction of LRO. 

Using the spin and charge distributions Q(i) — 1 (one 
electron per site) and S(i) = (—l) l Se z (AFM order in the 
arbitrary direction e z ), equation (^) yields two electronic 
bands characterized by the dispersion relations 



E 



i/v 



± 



(US) 2 



k e (-ir /2a, ir /2a] (19) 
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where = — 2t cos (ka) is the dispersion relation of non- 
interacting electrons and a is the lattice constant. Each 
of these levels is double-degenerate (a = ±1). Given the 
symmetry of the spectrum and the fact that only half 
the states are occupied, one can easily see that in the 
Hartree-Fock ground state, all the states in the valence 
band (E^. < 0) are occupied, while all the states in the 
conduction band (E^. > 0) are empty. The two bands are 
separated by the usual Mott-Hubbard charge transfer gap 
opened at the Fermi surface (k = ±n/2a), of magnitude 
2US. 

Using the valence band wavefunctions in Eqs. (^J) and 
([?]), we obtainlij the self-consistent spin amplitude 



V4 



(usy 



(20) 



where N is the number of sites and the sum is performed 
over the Brillouin zone k G (— 7r/2a, 7r/2a]. This equation 
has three solutions. One is trivial (S = 0). For 5^0 the 
equation depends only on 5" 2 . Consequently, the mean- 
field ground state is doubly-degenerate: both +S and —S 
satisfy Eq. (^0|) and give rise to self-consistent ground- 
states which differ from each other only through the fact 
that all the spins are flipped from one ground-state to 
the other one. Strictly speaking, the HF solution gives 
an infinite number of degenerate ground states, because 
the direction e z is arbitrary (this is a direct consequence 
of the fact that the Hubbard model is rotationally invari- 
ant). However, since a real chain is embedded in a 3D 
crystal, crystal field interactions will lift the rotational 
degeneracy, and fix one particular direction of orientation 
(easy axis) for the spins (for instance, along the chains). 
Consequently one particular direction e z is favored, and 
the mean-field ground state is doubly degenerate. 

Since all the states of the valence band are occupied, 
the energy of the HF ground state is simply given by: 



E GS = 2j2E v k +NU[S^ 



(21) 



where S is given by the self-consistency condition (20) 



The magnitude of the ground state energy per site, in 
units of t, is plotted in Fig. (|2|) as a function of U /t (the 
full line). The following features can be observed: in the 
U/t — > limit (noninteracting electrons), the energy of 
the ground state has, indeed, the expected value Egs —* 
—ANt/n. In the strong interaction limit U/t—* oo, the 
energy of the ground state goes to zero as expected (since 
in this limit every site is single occupied and virtual hop- 
ping is suppressed). For comparison, the exact ground- 
state energy obtained from the Bethe Ansatz (Eq. (|lS|)) 
is also shown (dashed line). The asymptotic value of the 
HF energy is found to be given by Egs —2Nt 2 /U. It 
is well known that in this limit, the Hubbard model is 
equivalent to an AFM Heisenberg modelJi3 with a cou- 
pling constant J =ri4i 2 /U , and that its true ground state 
energy per site isfl E GS -> -NJln2 = -2.77Nt 2 /U. 



This suggests that the Hartree-Fock method is a good 
starting point, from which to incorporate fluctuation cor- 
rections which lower the energy. 



2. Configuration Interaction treatment of the undoped chain 

The HF Slater determinant of the undoped AFM 
ground state is invariant to translations by 2a (AFM or- 
der must be preserved) . While it is possible to include in 
the CI set of Slater determinants excited HF states of the 
undoped chain obtained by exciting electrons from the 
valence to the conduction band, it turns out that lower 
energy self-consistent HF configurations can be gener- 
ated by breaking the translational symmetry of the un- 
doped mean-field AFM background. This is facilitated by 
considering excited states of the AFM background which 
can accommodate charge carriers in localized states deep 
within the charge transfer gap rather than within the 
Mott-Hubbard bands. The lowest energy self-consistent 
excited state is the undoped (neutral) domain-wall, which 
describes tunneling from one mean-field ground-state to 
the other mean-field ground-state. Since the AFM order 
rotates by 7r when crossing the domain wall, we consider 
either one domain wall on an odd-site chain, or a pair of 
domain walls on an even-site chain, in order to impose 
cyclic boundary conditions. Fig.|| depicts a typical self- 
consistent configuration containing two neutral domain 
walls, one centered between sites 10 and 11, and one cen- 
tered between sites 30 and 31. The charge Q(i) = 1 ev- 
erywhere. Near the domain wall the self-consistent spin 
magnitudes S z (i) adjust such that each neutral domain 
wall carries a spin 1/2 (with a projection in the same 
direction as its core spins), suggesting that this excita- 
tion is a neutral fermion. This is confirmed from the 
electronic structure shown in FigJJ. 

We obtain self-consistent configurations containing two 
neutral domain-walls at all possible distances from each 
other, either having opposite orientations (i.e. total chain 
spin 0) or same orientations (total chain spin ±1). Since 
the AFM ground state has total spin and states with 
different total spin do not mix, we need only include in 
the set of CI Slater determinants states of total spin 
zero, i.e. those having the neutral domain walls "paired" 
(with opposite orientations). Since all the possible config- 
urations with two neutral domain- walls have very similar 
energies, they must all be included along with the AFM 
undoped ground state in making up the variational trial 
wavefunction. We must also include both AFM undoped 
mean- field ground-states in the CI set. This can be easily 
seen from Figj|, where half of the chain is in one AFM 
ground state, and the other half is in the other (flipped) 
ground state. Therefore, this state will have equal over- 
lap with both AFM ground states, although the AFM 
ground-states themselves are orthogonal to each other. 
For an TV-site chain, the CI set contains a total of 2+N(N- 
2)/2 Slater determinants, two being the undoped AFM 
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ground states, and the rest being the N-2 possible states 
with paired domain- walls at different distances from each 
other, each of which can be translated N/2 times along 
the chain. 

The total ground-state energy found with the HF 
(circles), CI (squares) and BA (diamonds) methods for 
chains of different lengths N for U/t = 5 and U/t = 50 
are shown in Fig.||. For all three methods the total en- 
ergy of the chain is proportional to the length of the 
chain. The addition of the configurations with a pair 
of neutral domain walls in the CI method improves the 
ground-state energy considerably. It is obvious, however, 
that as the length of the chain increases configurations 
with two, three, four and more pairs of neutral domain 
walls should be included in the CI set in order to arrive 
at a perfect agreement with the exact Bethe Ansatz so- 
lution. It is interesting to remark that even if only the 
configurations with one pair of neutral domain walls are 
included, the non-degenerate CI ground state is such that 
(S z (i)) = for any site i of the chain, although the anti- 
ferromagnetic correlations remain very strong. This is a 
consequence of the fact that the CI set of Slater determi- 
nants contains an equal number of states with the spin 
at the site i up and down, so in average each spin expec- 
tation value is vanishing. Thus, the CI wavefunction is 
much more successful in mimicking the properties of the 
exact BA ground-state. 

B. Charged solitons in the doped ground state: the 
spin-bag and the charged domain-wall 

1. Hartree-Fock results 

If we numerically solve the HF equations (|],|]|lj) f° r a 
TV-site chain with N — 1 electrons, we find three types of 
charged self-consistent solutions: the spin-bag (or spin- 
polaron) (Fig||), the charged domain- wall centered on 
site (COS domain wall) (Fig. [7]) and the charged domain- 
wall centered between sites (CBS domain wall) (Figj|). 
The spin-polaron is created by trapping the hole in a 
small ferromagnetic core, which only affects the LR AFM 
order locally. The domain walls are topological excita- 
tions, since the AFM order is rotated by it as one goes 
through the domain wall. Therefore, in order to impose 
cyclic boundary conditions, we must consider an odd- 
site chain (or we may take an even-site chain and add 
two holes, leading to the appearance of a domain wall- 
anti-domain wall pair). The localization length of the 
hole decreases as U/t increasing for all three excitations. 
The spin and charge at sites far from the distortion equal 
the undoped mean-field ground-state values. 

The electronic spectra corresponding to the configu- 
rations shown in Figs. ||,|7] and || are shown in Figs. || 
and [lO. First panel in Fig. || corresponds to the un- 
doped ground state of a N = 40 sites chain. As dis- 
cussed before, the electronic spectrum consists of two 



bands of N = 40 states each. The valence band is com- 
pletely filled, the conduction band is completely empty, 
and there is a large charge transfer gap between them. 
Adding one hole on the same N = 40 sites chain, and 
keeping the cyclic boundary conditions, leads to the ap- 
pearance of the spin-bag shown in Fig.[| Its electronic 
structure is shown in the right panel of Figjjj. There is 
a localized level (n — 1) well below the valence band, 
the valence band contains 38 states, there are 3 localized 
levels deep inside the Mott-Hubbard charge transfer gap, 
and finally the conduction band also has 38 levels. Since 
there are N — 1 = 39 electrons in the system, only the 
localized level below the valence band and the valence 
band states are occupied. Since the valence band is spin 
paired (having an even number of states), this means 
that the total spin of this excitation is 1/2, associated 
with the spin of the electron on the localized level. The 
fact that the spin-bag carries a 1/2-spin is also easy to 
deduce from Fig. |[ because of the small ferromagnetic 
core. Thus, we conclude that the spin-bag is a charged 
fermion. 

The charged domain wall electronic structures are 
shown in Fig. [l0|, with the CBS domain wall in the left 
panel, and the COS domain wall in the right panel. In 
this case, we study a chain with N — 1 = 40 electrons and 
TV = 41 sites, so that we can impose cyclic boundary con- 
ditions again. We can see that in both cases there are 4 
localized levels inside the Mott-Hubbard gap, N — 1 = 40 
occupied states in the valence band, and 38 states in the 
conduction band. The localized levels of the COS domain 
wall are degenerate. The degeneracy is lifted for the CBS 
domain wall, and the upper discrete level is pushed quite 
close to the lower edge of the conduction band. How- 
ever, it is still a localized level (this is easily checked by 
plotting its wavefunction). In both cases we have a fully 
spin-paired valence band, and therefore the total spin of 
these excitations is zero. Since they carry the charge of 
the hole, the domain walls are charged bosons, in analogy 
to the charged solitons of polyacetylene.EZrES 

In order to establish the relevance of these different 
charged spin configurations, we compare their excitation 
energies (defined with respect to the undoped (half-filled) 
ground state) as a function of U/t in Fig. ( |ll| ) . As we 
can see, the domain walls are the low-energy excitations 
for U/t < 6.5, while the spin-polarons become the low- 
energy excitations for U/t > 6.5. As U/t — > 0, the core 
size of the domain- walls diverges roughly like t/U. As a 
result, in this limit the COS and the CBS domain-walls 
are very extended objects which become indistinguish- 
able and degenerate. However, as U/t increases the core 
becomes more and more localized, and the CBS domain- 
wall becomes energetically favorable relative to the COS 
domain- wall. 

However, this static HFA does not take into account 
the lowering of energy of these excitations due to trans- 
lations along the chain. From the simple inspection of 
the spin distributions of the spin-bag and of the domain 
wall, we can easily deduce that while a domain wall can 
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move freely along the chain, the spin-bag is rather im- 
mobile. Moving the center of the domain wall by one 
site (by interchanging the hole with the spin at the right 
or at the left) necessitates only some rearrangement of 
the magnitude of the core spins, while their orientation 
is automatically correct. Consequently, the domain wall 
lowers its energy by an amount of order t through hop- 
ping along the chain. However, if a spin-bag moves only 
by one site, the translated spin must be flipped (which 
would require the the total spin of the spin-bag to like- 
wise flip). In order to conserve its spin, the spin-bag 
must tunnel to the second nearest neighbor. This is a 
second order process, and consequently the spin-bag low- 
ers its energy only by an amount of order t 2 /U through 
motion. As we demonstrate below, using the Configu- 
ration Interaction method, these qualitative arguments 
are valid. When soliton dynamics is incorporated, it is 
the charged bosonic domain wall that proves to be the 
relevant charged excitation of the Hubbard model for all 
values of U/t. 

2. Configuration Interaction treatment of the spin-polaron 

Consider a spin-polaron on a chain with 2N sites (the 
number of sites is even so that we can impose cyclic 
boundary conditions). Using the CI method we evalu- 
ate the kinetic energy of the mobile, charged spin-bag. 
As suggested above, we only need to include in the set of 
Slater determinants configurations translated by an 
even number of sites from the initial HF configuration. 
Let l^e-uen), \^ odd) be the HF determinants for the spin- 
polaron centered at an even and odd site, respectively, 
and let S z = £\ = § J2i, a ac \a c «y be the total spin 
operator in the ^-direction. Then, S z \^! even ) = ^|\& even) 
while S z \^ oc id) = — \\^odd) (or viceversa), and therefore 

odd\^ even) — 0. Since the Hubbard Hamiltonian com- 
mutes with S z , it follows that dd\H\^ even) — 0. From 
the CI equation ([n]) we conclude that there is no mixing 
between states with the spin-polaron on one sublattice 
and states with the (opposite spin) spin-polaron on the 
other sublattice. Therefore, on a chain with 2N sites we 
only need to mix N Slater determinants in order to obtain 
the spin-polaron ground-state within the CI method. 

If the initial self-consistent HF spin-polaron config- 
uration l^i) is composed of the one-particle occupied 

orbitals (f>p(i,a), the one-particle orbitals of the state 
|\& n +i) translated by 2na will simply be chosen as 

</>p™ +1 ^ (i,o-) — <f>p (i — 2n, a) (cyclic boundary conditions 
are assumed). The overlap matrices O nm and Tt nm are 
then calculated and the CI matrix equation ([fl]) solved. 
Numerically, the largest amount of time is spent com- 
puting the Tinm matrix elements. Due to various sym- 
metries, there are only N/2 distinct matrix elements. 

Given the structure of the CI equation (|ll|), we can 
readily see that its solutions are of the form 



N 

l**) = X> (ttna, °l*») ( 22 ) 

n=l 

where a is the lattice constant and there are N distinct 
k values. These values satisfy the periodicity condition 
exp (ikN2a) = I since translating any spin-polaron con- 
figuration by the total chain length 2Na leaves the config- 
uration unchanged. Therefore, the distinct wavevectors 
are k = mir/Na, m = 0, 1, N — 1. The reduced Bril- 
louin zone [0, 7r/a) (or, symmetrically, the (— ir/2a, ir/2a\ 
interval) is due to the motion of the spin-polaron on 
only one sublattice, and corresponds to states with spin- 
up (for instance). The band corresponding to the spin- 
polaron on the other sublattice will have the same struc- 
ture, but corresponds to spin-down states. In the end 
we recover the typical spin-degenerate band expected for 
fermions. 

Given the general form of the wavefunction, the dis- 
persion relation of the spin-polaron follows from the ex- 
pression 

_ (*E fc |K|tt fc ) _ £? =1 exp(2iA: a (n-l))ft ln 
(*fc|**> J2" =1 exp(2ika(n-l))0 ln 

In deriving the last equation, we used the symmetry 
properties of the matrices TL nm an d O nm , namely that 
the (nm) matrix element only depends onn-m. Strictly 
speaking, E(k) is the energy of the whole chain con- 
taining the spin-polaron and will strongly depend on the 
length of the chain. We extract the dispersion relation of 
the spin-polaron from a fit of the form 

E{k) = 2Ne GS + E pol {k) (24) 

Here ecs is interpreted as the ground-state energy per 
site of the undoped AFM (for a very long chain, most of 
the sites are unaffected by the presence of the single spin- 
polaron). We define E po i{k) as the dispersion relation of 
the spin-polaron itself. In other words, the energy of 
the spin-polaron is defined as the difference between the 
energy of the chain with the spin-polaron, and that of an 
undoped chain. 

We plot E po i(k) versus k in Fig. O, for U/t — 5 and 
chains of various lengths. The various curves fall on top 
of each other, thus proving that the fit ( p4| ) is legitimate. 
Also shown is the excitation energy of the static hole- 
doped spin-polaron E^ f (the full line) , as obtained from 
the Unrestricted HFA (also defined with respect to the 
energy of the undoped chain) . Clearly, translation lowers 
the total energy of the spin-polaron, with the most stable 
state corresponding to k = ix/2a. The total kinetic en- 
ergy gained is, however, only of the order of t 2 /U . This 
is shown in Fig. O, where we plot both the kinetic en- 
ergy gained E po i[^) — E^J/ (circles) and the width of 
the spin-polaron band, E po i{^-) — E po i(0) (squares), as a 
function oit 2 /U. The log- log graph is linear with a slope 
of unity as expected, since the spin-polaron must tun- 
nel two sites (second order hopping process) to the next 
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spin-allowed position. Clearly, this charged fermionic ex- 
citation is relatively immobile. 

We conclude that in the large U/t limit the CI cor- 
rection to the spin-bag energy is negligible, due to the 
immobility of this excitation. As a result, the energy 
of the spin-bag varies with U/t as shown in Figjll], for 
large U/t, and it saturates above — 1.5t as U/t — » oo. 
As already discussed, it is known that in the U/t — ► oo 
limit the energy of the doping hole is — It. This discrep- 
ancy suggests that the spin-bag does not provide a good 
description for the charge carrier. 

3. Configuration Interaction treatment of the charged 
domain-wall 

To investigate an isolated charged domain wall, we 
consider chains with an odd number (2N + 1) of sites. 
As shown in Figs, [t] and ||, there are two types of self- 
consistent charged domain-walls, namely the COS (cen- 
tered on site) domain wall and the CBS (centered be- 
tween sites) domain wall. Since the charged domain walls 
are bosons (S z \^>) = 0), there is non- vanishing overlap 
between states with the domain- wall centered on different 
sublattices. Unlike the fermionic spin-polarons, we must 
include all possible translations in the CI Slater determi- 
nant set I*,-). For a 2iV+l-site chain, there are 2(2iV+l) 
distinct COS domain wall configurations, and 2(2A" + 1) 
distinct CBS domain wall configurations. The reason for 
the factor 2 is that translation of a domain wall by 2N + 1 
sites takes it into a domain wall centered at the same site 
as in the initial configuration, but with all spins flipped, 
due to the 7r-difference in the AFM ordering of the spins 
on the two sides of the domain wall. A second translation 
around the whole chain is necessary in order to regain 
the initial configuration. As a result, there is a four-fold 
increase in the number of possible configurations for a 
domain wall, as compared to a spin-polaron on a chain 
of almost the same length. We generate the translations 
with an even number 2n of lattice sites in the same way 
as for the spin-polaron, cf>^ n+1 ' (i, a) — <j>p {i — 2n, a) 
if < n < N (first translation around the chain) and 

4 2 " +1) («,ct) = <Pp\i ~ 2n,-a) if N < n < 2N (sec- 
ond translation around the chain). Here, we remind the 

reader that <fip n \i,a) represents a particular (p) occu- 
pied one-electron orbital of the static Hartree-Fock Slater 
determinant \^ m ) which places a static magnetic soli- 
ton at site m. For translations with an odd number 
2n — 1 of sites, the wavefunction mapping is 4>p 2n ^ (i, a) — 

<f)p (i — (2n — I), — a) if I < n < N (first translation 

around the chain) and (f> p (i, a) = (f>p (i — (2n — 1), a) 
if N < n < 2N+1 (second translation around the chain). 

Let us first consider only including one type of domain 
wall (either COS, or CBS), in the CI wavefunction. In 
this case, we can again conclude that the solutions of Eq. 
( pr| ) must be of the form 



2(27V+1) 

|**)= E e' fe "°l*n) (25) 

where \^ n ) is the configuration translated by n — 1 
sites from the initial HF configuration The pe- 

riodicity condition is now e tk2 ( 2N + 1 ) a — i an d the al- 
lowed values of k are given by k = mir/(2N + l)a,m = 
0, 1, ...,2(2iV— 1) — 1. Clearly, the domain- wall dispersion 
band is extended over the full Brillouin zone [0, 2ir) (or 
the symmetric version (— tt,tv]). The dispersion relation 
is given by 

E (k) = <**M* fc > = Sn^'^PW"-^ 

(*k\* k ) E' ( =f +1) exp^(«-l)a)O nl 

(26) 

As in the case of the spin-polaron, we extract the disper- 
sion relation of the domain-wall from E(k) by subtracting 
the energy of the undoped chain 

E dw (fc) =E(k)- (2N + l)e GS (27) 

In Fig. [l4| we show the dispersion relations Ed w (k) 
vs. k for both CBS (left panel) and COS (right panel) 
domain walls on chains of different length 2N + 1, and 
U/t = 5. The excitation energy of the static configura- 
tion (obtained from the unrestricted HF search) is also 
shown. Again, various dispersion curves fall on top of 
each other, validating the fit of Eq. (p7[). Comparing Fig. 
[l4| with Fig. |lj, it is immediately apparent that the dis- 
persion band of the domain walls is much wider. In fact, 
the band of the COS domain wall extends up to 4t (not 
shown). Comparing the bottom of the dispersion band 
to the static HF excitation energy of domain wall (shown 
as a full line), we see that translational motion lowers the 
energy of the domain wall by about t (as opposed to only 
0.3i when U/t — 5 for a spin-polaron). While the bot- 
tom of the dispersion band is basically identical for both 
types of domain walls, the top is very different. Excited 
states with energy E{k) > require the incorporation of 
the excited state configurations of the single Slater deter- 
minant (from static Hartree-Fock) in the CI set. Clearly 
these excited state, static configurations have energies 
comparable to the moving domain wall at high energy 
parts of the dispersion curve. The most likely candidates 
are those configurations in which electrons from the top 
of the valence band are excited into the bound discrete 
levels. If only one such excitation takes place, the en- 
ergy of the static configuration is raised by « U/2 (the 
difference between the energy of the level at the top of 
the valence band, and that of the first empty localized 
level). By mixing such configurations in the CI Slater 
determinant set, we obtain modifications to the upper 
part of the dispersion relation, while the bottom remains 
unchanged. Since we are interested in the kinetic energy 
gained by the domain wall through translation {E(k) < 
region), we will neglect these higher-energy processes in 
what follows. 
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A technical issue that emerges is the effect of mix- 
ing both the COS and the CBS domain wall configura- 
tions when calculating the CI wavefunction. While one 
might hope for an improvement in the overall energy for 
the mobile, charged soliton, this is not the case. The 
reason is that each set of configurations by itself gen- 
erates basically the same CI wavefunctions \^k) at the 
bottom of the dispersion band rather than linearly in- 
dependent ones. This can easily be seen numerically if 
we analyze the eigenvalues of the overlap matrix O nm . 
Suppose that A is an eigenvalue of this matrix, and that 
(ofi)i=i,AT is the corresponding eigenvector (for simplicity, 
we use N as the dimension of the overlap matrix) . Defin- 
ing \*f?) = J2i=i a i\^i) it is straightforward to show that 
(^1^) = A. When we mix both sets of configurations to- 
gether, we find many vanishing eigenvalues A = 0, which 
imply I*) = (numerically, we use the Singular Value 
Decomposition technique as a diagnostic for vanishing 
eigenvalues). This proves that there are redundant lin- 
early dependent combinations in the set of Slater deter- 
minants l^i). We remove these linearly dependent states 
to find the CI ground-state. In particular, by mixing 
COS and CBS domain wall configurations, the resulting 
low-energy spectrum is the same as that found by using 
only the lower energy CBS domain wall configurations. 

The previous analysis gives us the lowest energy of a 
single hole (charged domain wall soliton) on the chain 
that can be obtained within the CI approximation. In 
comparing this energy to the one obtained from the ex- 
act Bethe Ansatz solution, it is not appropriate to di- 
rectly compare the total chain energies. The reason, as 
already proved, is that there is a large contribution to 
these energies proportional to the number of sites in the 
chain, the proportionality constant being the undoped 
ground-state energy per site (see Eq. (|27|)). The HFA 
gives a higher undoped energy per site than the exact 
BA energy (see Fig. ||), and the CI approach does not 
improve it unless we also add states with pairs of un- 
charged domain- walls. Our aim is to isolate the energy 
of the doping hole. Therefore, we compare Ed w {^) (the 
lowest CI energy of the domain wall itself) with the cor- 
responding doping hole energy extracted from the Bethe 
Ansatz. This comparison is shown in Fig. [l5|. In order to 
find the doping hole energy from the Bethe Ansatz, we 
evaluate the exact ground-state energy of a chain with 
2N + 1 sites and 2N electrons ( half M = N of which 
have spin-down), for various values of N. This set of 
energies is seen to be well fitted by an expression of the 
form E{N) = (2N + l)ecs + E 0l where ecs is in excel- 
lent agreement with the BA ground-state value predicted 
by (|l8|). As in the CI approach, we identify Eq with the 
energy of the hole. The BA energies of the hole, as a 
function of U/t are shown as squares in Fig. ^|. They 
indeed go to — 2t in the U/t — > oo limit, as expected. In 
the U/t — > limit, the energy of the doping hole is ex- 
pected to go to zero, since in this U = limit the system 
is a metal. 



The CI domain- wall energies Ed w = Ed w ( — ) are shown 
as full circles in Figure [l5|. The agreement with the Bethe 
Ansatz energy is striking. For U/t > 5, the fit ( pjj ) is 
excellent and the error bars on the domain wall energies 
are extremely small. However, as U/t — > 0, the size of 
the domain wall increases significantly (it is around 20 
sites for U/t = 2) and therefore extremely long chains 
need to be considered for a good fit. The two upper lines 
correspond to the static HF energies obtained for a self- 
consistent CBS domain wall (triangle down) and a COS 
domain wall (triangle up). The diamonds show the CI 
results for the spin-polaron. Clearly, the translational 
motion of the domain wall (included in the CI approach) 
drastically lowers its overall energy. The kinetic energy 
saved is of the order t over most of the U/t parameter 
range. 

The agreement between the domain- wall energy as cal- 
culated in the CI approach and the exact doping-hole 
energy as given by the Bethe Ansatz is quite remark- 
able, over the whole range of U/t parameters. The only 
disagreement appears for U/t < 2, where the domain- 
walls become extremely delocalized and the numerical 
calculations are very difficult. The CI solution is not ex- 
act because the HF description neglects the presence of 
additional neutral domain-walls pairs in the AFM back- 
ground. While the addition of such pairs improves the 
accuracy of the CI method relative to the exact solu- 
tion, it makes the calculation more cumbersome, due to 
the large increase in the number of possible configura- 
tions. When the contribution of this background is re- 
moved, the CI solution of a single charged domain wall 
moving around the chain agrees very well with the exact 
solution. This suggests that the renormalization of the 
charged domain wall energy by neutral domain wall pairs 
is relatively small. The agreement with the Bethe Ansatz 
is not limited to the bottom of the charged domain wall 
dispersion curve. In the large U/t limit, the domain-wall 
dispersion band is indeed given by Ed w (k) — 2tcos(ka), 
as requiredJiil For example, the dispersion relation for 
domain-wall corresponding to U/t — 100 is shown in Fig. 
|l^ . At such high U/t values, the typical energy of the 
configurations containing electrons excited on the mid- 
gap levels is of order {7/2. They do not influence the 
lowest domain wall band. As U/t is decreased, these ex- 
cited configurations simply modify the high energy part 
of the domain- wall dispersion relation. 

In conclusion, the mobile charged bosonic domain wall 
excitation is the relevant charged excitation of the ID 
chain described by the Hubbard model for all values U/t. 
Although the static charged spin-polaron has lower exci- 
tation energy than the static domain- wall for U/t > 6.5, 
when quantum dynamics is taken into consideration (CI 
method) the mobile charged domain wall turns out to be 
the lowest energy excitation. 
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C. Spin-charge separation 

A particularly striking effect can be recaptured in the 
CI method if one adds a hole to an even-site chain. In the 
HFA, this leads to the appearance of one spin-polaron, 
since a single charged domain wall is incompatible with 
the cyclic boundary conditions involving an even num- 
ber of sites. In the CI method, this charged spin-polaron 
is unstable to dissociation into a pair of more mobile 
domain- wall solitons. In particular the charged spin-bag 
dissociates into a charged bosonic domain- wall and a neu- 
tral fermionic domain wall. The translational kinetic en- 
ergy saved by the domain wall motion more than offsets 
the additional exchange energy cost in creating a pair of 
solitons from a single spin-polaron. 

We demonstrate this spin-charge separation effect us- 
ing the set of all the possible configurations with a pair of 
a charged and a neutral domain walls in the CI method. 
Again, since states with different spins do not mix, we 
need only keep configurations in which the uncharged 
domain wall has the same spin, either +1/2, or -1/2. As 
in the undoped case, configurations with a pair of domain 
walls connect the two possible AFM ground states. This 
leads to a considerable change in the background energy 
(see Fig.||). As a result, we define the energy of the dop- 
ing hole, in this case, with respect to the CI energy of an 
undoped chain with a pair of neutral domain-walls. This 
allows us to properly account for the renormalization of 
the background energy of the chain by the pair of do- 
main walls. Physically, we interpret this in the following 
way. The true ground state of the undoped chain has a 
certain number of pairs of neutral domain walls. When 
the chain is doped, the doping hole is bound into one of 
the already existing neutral domain walls, transforming 
it into a charged bosonic domain wall. All the other pairs 
of neutral domain walls remain largely unaffected. 

The results of this analysis are plotted in Fig.|l^. We 
use the Bethe equations to calculate the ground-state en- 
ergy of even-site chains with one doping hole. The BA 
excitation energy for the doping hole added to an even- 
site chain is shown as circles. The exact chain energy is 
found to be well fitted by E(N) = Ne GS + E , with N 
being the even number of sites of the chain, and ecs the 
undoped ground-state energy per site. Since Necs is the 
total energy of the undoped chain, we again identify E Q 
as the energy of the doping hole. The energy of a charged 
domain wall - neutral domain wall pair, obtained using 
the CI method for even-site chains, is shown as squares. 
Clearly, the CI method recaptures the physics as well 
as the energetics of the BA solution to a high degree of 
accuracy. Even closer agreement can be achieved if con- 
figurations with more pairs of uncharged domain-walls 
are included in the CI basis set. This result provides a 
very clear illustration of the spin-charge separation phe- 
nomenon known to exist in the ID Hubbard model. Upon 
doping, the relevant charge excitation is not the quasi- 
particle excitation (spin-polaron) which carries both the 



spin and the charge of the doping hole, but rather a de- 
confined pair excitation consisting of a charged bosonic 
domain-wall (carrying charge but no spin) and a neutral 
fermionic domain- wall (carrying spin but no charge). 

Another interesting result of this analysis is that the 
excitation energy needed to add a hole to an even-site 
chain is different from that needed to add a hole to an 
odd-site chain (see Fig.|l7], circles and diamonds). Such a 
difference would be irreconcilable if the excitation was a 
local quasiparticle, in which case different boundary con- 
ditions are expected to introduce variations of the order 
0(1/N). However, in the presence of the charge-spin sep- 
aration, there is a simple interpretation. A doped chain 
with an odd N number of sites has an even N — 1 number 
of electrons. Its ground state has zero spin (all electrons 
are paired). This is well described by a single charged 
domain-wall. On the other hand, the doped chain with 
an even N number of sites has an odd N — 1 number 
of electrons. In order to represent the total spin 1/2 of 
the unpaired electron, it is necessary to include both a 
charged domain wall and a neutral spin- 1/2 soliton. We 
associate the difference AE between the excitation en- 
ergy of a single hole on odd and even-site chains with the 
energy of the additional neutral spin- 1/2 soliton. This 
interpretation may be verified independently by evaluat- 
ing the energy of a single neutral spin-1/2 domain wall 
on an odd-site, undoped chain. The BA ground-state en- 
ergy of the undoped odd-site chain varies with the odd- 
N number of sites as E(N) = Necs + Ed w , where eas, 
the undoped ground-state energy per site, has the same 
value as obtained from fits of even-site chains. We iden- 
tify Edw as the energy of the neutral spin soliton. The 
energies Ed w and AE are compared in Fig.|l8|. The good 
agreement between these independent measures of energy 
confirms that the spin soliton is a well defined concept 
even in the presence of doping. The existence of the spin 
soliton accounts for the difference in energy between odd 
and even site chains even in the thermodynamic limit. 
It facilitates our identification of the spin soliton in the 
Bethe ansatz with the neutral spin-1/2 domain wall of 
the CI approximation. 

V. CONCLUSIONS 

In this paper we have demonstrated the validity of the 
Configuration Interaction Method as a technique for go- 
ing beyond mean field theory in strongly correlated elec- 
tron systems. The method recaptures the essential phys- 
ical features as well as energetics of the exact solution of 
the ID Hubbard chain. The CI method provides a sys- 
tematic way to improve and go beyond the Hartree-Fock 
Approximation, by incorporating essential quantum dy- 
namics and tunneling effects which are excluded from the 
mean-field theory. 

We showed that a charged bosonic domain wall can 
lower its kinetic energy by about t for all U/t values, 
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while the immobile charged spin-polaron can only lower 
its kinetic energy by an energy of the order of t 2 /U. As 
a result, the mobile charged bosonic domain wall is the 
low energy charged excitation of the Hubbard chain for 
all values of U/t, and its excitation energy and disper- 
sion band are in good agreement with the predictions of 
the exact Bethe-Ansatz solution. We also showed that 
it is energetically favorable for the (quasiparticle-like) 
charged fermionic spin-bag to dissociate into a charged 
bosonic domain wall (which carries the charge but no 
spin) and a neutral spin- 1/2 domain wall (which carries 
the spin but no charge). Clearly the CI method recap- 
tures the physics of spin-charge separation known to exist 
in the ID Hubbard model. It may be fruitful to perform 
a more detailed study of the interaction between various 
types of domain walls, both charged and uncharged. Var- 
ious correlation functions and response functions can also 
be calculated for detailed comparison with exact results. 

Given the validity and accuracy of the Configuration 
Interaction Method for the ID problem, we believe that it 
provides a natural route to describe effects beyond mean- 
field theory in the doped ground-state of the 2D Hubbard 
model and other models of strongly correlated electrons 
in higher dimensions. We present this study elsewhere.u 
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FIG. 4. Electronic structure of the 40-site chain with the 
two neutral domain walls shown in Fig.[| Each domain wall 
has 4 discrete levels bound in its core. The spins on the two 
occupied bound levels are oriented in the same direction as the 
core spins of the domain wall. In the configuration shown in 
Fig.^ , there is a +1/2 and a -1/2 domain wall, and therefore 
all levels are spin paired and the total spin of the chain is zero. 
However, for two +1/2 (-1/2) domain walls, all the occupied 
discrete levels have +1/2 (-1/2) spins, and the total spin of 
the chain is +1 (-1). 

FIG. 5. Ground-state energy (in units of t) of a chain of 
size N, calculated with the HF (circles), CI with two neutral 
domain walls (squares) and BA (diamonds). The left panel 
corresponds to U/t = 5, while the right one corresponds to 
U/t — 50. Although the energy scale is very different, in 
both cases the CI method significantly improves the agree- 
ment with the exact Bethe Ansatz solution. 
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FIG. 6. Self-consistent charge (upper line) and spin distri- 
butions for a charged spin-bag on a 40-site chain, for U/t = 5. 
The spin-bag is a charged fermion. 

FIG. 7. Self-consistent charge (upper line) and spin distri- 
butions for a domain wall centered on site (COS) on a 41-site 
chain, for U/t = 5. 

FIG. 8. Self-consistent charge (upper line) and spin dis- 
tributions for a domain wall centered between sites (CBS) on 
a 41-site chain, for U/t — 5. The domain walls are charged 
bosons. 

FIG. 9. Electronic spectra for an undoped AFM chain 
with 40 sites (left panel), and for a 40-site chain with a 
charged spin-bag (right panel). U/t = 5. 



FIG. 10. Electronic spectra for a chain 41 sites chain 
with a CBS domain-wall (left panel), and a COS domain- wall 
(right panel). U/t — 5. 



FIG. 14. The dispersion relations Ed w (k) vs. k for both 
CBS (left panel) and COS (right panel) domain walls on 
chains of different length 2N + 1 = 17, ...,23, and U/t = 5. 
The excitation energy of the static configuration (obtained 
from the unrestricted HF search) is also shown as a full line. 
The extra kinetic energy gained through translation by the 
domain-wall is of the order of t. 

FIG. 15. Excitation energy, in units of t, for a mo- 
bile charged domain wall (circles) and a mobile charged 
spin-polaron (diamonds), as obtained from the CI approach. 
The exact excitation energy given by the Bethe-Ansatz 
method is shown by squares. The domain-wall CI energy is in 
excellent agreement with the exact BA results (also see inset), 
while the spin-polaron CI energy is significantly different. For 
comparison, we also show the excitation energies for the COS 
and CBS domain walls as obtained from the static HFA (up 
and down triangles), proving again that the extra kinetic en- 
ergy gained by the moving domain wall is of order t for most 
U/t values. In contrast, the extra kinetic energy gained by 
the spin-polaron is of order t 2 /U — > as U/t increases, so in 
the large U/t limit there is almost no difference between the 
HF and CI results for the charged spin-polaron. We conclude 
that the charged domain-wall is the relevant excitation for all 
values of U/t. 



FIG. 11. Excitation energy in units of t of the charged 
spin-polaron, CBS and COS domain walls, as a function of 
U/t. The excitation energy is defined with respect to the un- 
doped chain. Within the HFA, for U/t < 6.5 the charged do- 
main-walls are the low-energy charged excitations, while for 
U/t > 6.5 the spin-polarons are the low-energy excitations. 
However, the HFA approximation neglects the kinetic energy 
gained by these charged excitations through translation along 
the chain. When this is taken into account within the Config- 
uration Interaction approximation, the mobile domain-wall is 
found to be the low-energy charged excitation for all values 
of U/t (see Fig. 

FIG. 12. Dispersion band for the spin-polaron, E po i(k) 
vs. k, with Ep i(k) extracted from Eq. (E4) for chains of 
length 2JV = 14, 16, 22 sites and U/t = 5. Also shown is 
the excitation energy of the static hole-doped spin-polaron 
Epof (t ne full line), as obtained from the Unrestricted HFA. 
Translation lowers the total energy of the spin-polaron, with 
the most stable state corresponding to k = n/2a. However, 
the kinetic energy gained through translation is quite small. 

FIG. 13. The kinetic energy gained by the delocalized 
spin-polaron E po i(j^) — E po f (circles) and the width of the 
spin-polaron band, E po i(j^) — E po i(0) (squares) as a function 
of t 2 /U. The log- log graph is linear with a slope of unity 
as expected. The spin-polaron must tunnel two sites to the 
next allowed position, through a second order hopping pro- 
cess. This charged fermion is rather immobile. 



FIG. 16. The dispersion relation for a domain-wall on a 
chain of length 2iV+l = 17, 19, 21 and U/t = 100. In the large 
U/t limit the disperaien relation of one single hole is given by 
E dw {k) = 2tcos(ka)JB This is indeed in very good agreement 
with the dispersion band of the domain wall, proving again 
that this is the relevant charged excitation of the Hubbard 
chain. 

FIG. 17. Excitation energy (in units of i) for a charged 
domain- wall - neutral domain- wall pair, as obtained from the 
CI approach for even-site chains (squares). The energy of the 
doping hole added to an even-site chain, as obtained from 
the exact Bethe Ansatz solution, is shown by circles. Again, 
there is good agreement between the two methods. For com- 
parison, we also show the energy of the doping hole added to 
an odd-site chain, as obtained from the Bethe Ansatz (dia- 
monds) and the energy of an isolated charged domain wall on 
odd-site chains (triangles). These last two sets are the same 
as in Fig.|l5[ 

FIG. 18. The excitation energy of a neutral domain-wall, 
as obtained from the Bethe ansatz for undoped odd-site chains 
(circles). This agrees with (diamonds) the difference in the 
excitation energy of a single hole on an even-site vs. odd-site 
chain (shown individually as circles and diamonds of Fig.^) . 
The dotted line is the asymptotic fit Qt 2 /U. 
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